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ABSTRACT 

We consider a universe in which inflation commences because of a pos- 
itive cosmological constant, the effect of which is progressively screened by 
the interaction between virtual gravitons that become trapped in the ex- 
pansion of spacetime. Perturbative calculations have shown that screening 
becomes non-perturbatively large at late times. In this paper we consider 
effective fleld equations which can be evolved numerically to provide a non- 
perturbative description of the process. The induced stress tensor is that of 
an effective scalar field which is a non-local functional of the metric. We use 
the known perturbative result, constrained by general principles and guided 
by a physical description of the screening mechanism, to formulate a class 
of ansatze for this functional. A scheme is given for numerically evolving 
the field equations which result from a simple ansatz, from the beginning 
of infiation past the time when it ends. We find that infiation comes to a 
sudden end, producing a system whose equation of state rapidly approaches 
that of radiation. Explicit numerical results are presented. 
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1 Introduction 



Perturbation theory is an immensely gratifying tool. It almost always pro- 
vides quantitative answers for how a known system changes with the inclu- 
sion of a small, new effect. The great frustration in using the technique is 
that its answers become unreliable precisely when they are most interesting: 
when the new effect causes major changes. Our recent study of the quan- 
tum gravitational back-reaction on inflation illustrates both the utility of 
perturbation theory, and the frustration of not being able to push it further. 

The unperturbed system in our case is classical general relativity, the 
Lagrangian for which is: 

Here G is Newton's constant and A is the cosmological constant, assumed 
positive. On a spatially flat manifold the invariant element for a homogeneous 
and isotropic universe can be written in co-moving coordinates: 

g^,^,{t, x)dx''dx'' = -df + e'^^^^^dx ■ dx . (2) 

And the classical solution is: 



6dass(t) = Ht , (3) 

where H = (A/3)^/^ is the Hubble constant. If we specialize to the manifold 
X 3?, where each of the coordinate radii is H~^, then the 3- volume: 

V{t) = iZ-^e^^* , (4) 

is finite but grows exponentially. 

The perturbation we seek to study is the gravitational interaction between 
virtual infrared gravitons that become trapped in the expansion of spacetime 
and get pulled apart. Although we have computed this exactly at the lowest 
non-trivial order in perturbation theory |T|, 0], an intuitive understanding 
of the effect is necessary if we are to abstract it beyond the perturbative 
regime. The physical picture is that virtual gravitons of sufficiently long 
physical wave length are torn apart by inflation. This is the phenomenon 
of superadiabatic amplification, first studied by Grishchuk in 1974 Al- 
though infrared gravitons are continually produced in this way, the volume 
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of space expands so rapidly that the energy density of these gravitons re- 
mains a constant — and rather small — fraction of A/(87rG).[] However, as 
each graviton pair recedes, the intervening space is filled by their long range 
gravitational potentials. These potentials persist even after the gravitons 
that engendered them have reached cosmological separations. As new pairs 
are ripped apart, their potentials add to those already present. This is a 
secular effect and it obviously continues as long as inflation does. Because 
gravity is attractive the effect tends to counteract inflation — and hence to 
screen the cosmological constant. 

One might expect similar results from quanta other than gravitons but 
this is not so. To experience superadiabatic amplification a particle must 
be effectively massless with respect to H, and it must not possess classical 
conformal invariance [0, 0.0 One or the other of these two conditions ex- 
cludes every other known particle and most of the conjectured ones. The 
only contender, besides the graviton, is a massless, minimally coupled scalar. 
These do experience superadiabatic amplification, but global conformal in- 
variance prevents them from inducing a gravitational interaction comparable 
to that of gravitons f^. One might get a strong effect if such a scalar had 
non-derivative self-interactions, but it is difficult to understand why these 
would not also induce a substantial mass. 

Screening affords a simple and satisfying reformulation of inflationary 
cosmology and a beautiful resolution to the associated problems of fine tun- 
ing. Inflation starts, in this scheme, because the cosmological constant is 
positive and not unreasonably small. Inflation eventually ends due to the 
self-gravitation of virtual gravitons which have become trapped in the super- 
luminal expansion of spacetime. Inflation lasts for a very long time because 
gravitational interactions are weak, even at the GUT scale. One can be indif- 
ferent about adding matter because gravitons are the unique phenomenolog- 
ically viable quanta which induce screening. The only thing to avoid, in this 
scheme, is introducing an inflaton field and fine tuning its potential! Best 
of all, the infrared character of the screening mechanism means that it can 

^It is easy to show that there is on average one infrared graviton per Hubble volume. 
Even for inflation on the GUT scale this is only about 10^^^ of the energy density of the 
cosmological constant. 

^Massive particles are short range, so their virtual quanta seldom get far enough apart 
to become trapped in the expansion of spacetime. Conformally invariant particles are 
incapable, locally, of distinguishing between the conformally flat classical background (||) 
and flat space. 
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be studied reliably using quantum general relativity, in spite of the theory's 
lack of perturbative renormalizability and without regard to what happens 
at the Planck scale .0 

But there is a problem: the quantum gravitational back-reaction can only 
be studied perturbatively so long as it is weak. This regime is not without 
interest. For example, perturbative analysis shows that inflation lasts a long 
time and that conventional matter is incapable of competing with the 
quantum gravitational back-reaction 0. If one assumes a sudden end to 
inflation — which is certainly supported by the perturbative results — 
then it should be possible to predict the spectrum of density fluctuations in 
the perturbative regime. However, the most interesting questions lie frus- 
tratingly beyond the point where perturbation theory is valid. Hence the 
need for a non-perturbative model. 

Our strategy for creating such a model is to infer the induced stress tensor: 

T^,A9\ = ^ (Rf^u - ^9f.uR + 9t,u^ , (5) 

as a non-local functional of the metric which correctly reproduces the known 
perturbative effect and which captures the physical origin of screening gener- 
ally. Of course there is some ambiguity in this, but surprisingly little of any 
significance. Given an ansatz we can numerically integrate the field equa- 
tions as far into the future as is desired. What we find, for a simple ansatz, is 
that inflation ends suddenly over a period of about five e-foldings, following 
which the equation of state rapidly approaches that of pure radiation. 

This paper consists of eight sections, of which the first is drawing to a 
close. In Section 2 we show that, for the purposes of cosmology, the induced 
stress tensor can be parameterized as that of an effective scalar field which 
is a non-local functional of the metric. We also derive what this functional 
must be during the perturbative regime. In Section 3 we enumerate six 
principles which constrain the scalar functional generally. Section 4 gives a 
semi-quantitative model of screening which is of course the ultimate physical 
motivation for the choice of scalar. We discuss various ansatze for the scalar 

■^Infrared phenomena can always be studied using the low energy effective theory. This 
is why Bloch and Nordsieck were able to resolve the infrared problem in QED before 
the theory's renormalizability was suspected. It is also why Weinberg (t) was able to 
give a similar resolution for the infrared problem of quantum general relativity with zero 
cosmological constant. And it is why Feinberg and Sucher [|| were able to compute the 
long range force induced by neutrino exchange using Fermi Theory. 
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in Section 5. Section 6 describes our scheme for numerically integrating the 
dynamical system resulting from a simple ansatz. We also report explicit 
results. In Section 7 we reconstruct the potential of the effective scalar, 
analytically for large or small values of the scalar and numerically for any 
value. Section 8 is a discussion of our results. 

2 Effective scalar stress tensor 

The point of this section is to show that, for the purposes of cosmology, we 
can model the induced stress tensor (^) as that of a scalar field which is itself 
a non-local functional of the metric: 

T,A9] = d,<P[g]dM - 9,. {^g'^'dpMdM + • (6) 

We will also show that, since the potential can always be chosen to enforce 
conservation, one really needs only the functional (l)[g\. And we will use the 
known perturbative results to derive what cplg] must be when specialized 
to the classical background (H). 

There are three senses in which one might discuss the induced stress 
tensor, or any other functional of the metric. The first sense is generally for 
an arbitrary metric; the second is as a functional of h{t) for a spatially fiat, 
homogeneous and isotropic metric; and the third is as an explicit function 
of time for the case of perturbative corrections. When the basic symbol 
appears unadorned we mean the general quantity; the presence of a hat 
implies specialization to the spatially flat, homogeneous and isotropic form; 
and the perturbative version of the same quantity is denoted by a tilde. For 
example, a general invariant element is represented thus: 

ds^ = g^^{t,x)dx^dx'' . (7) 
Specializing to fiat, homogeneous and isotropic spacetimes gives: 

ds^ = -df + exp[2b{t)]dx ■ dx . (8) 
And we express the perturbative result as follows [0: 

ds^ = -dt^ + + A{t) exp [2Ht {l + D{t)}] , (9) 
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where perturbative expansions for the functions D{t) and A{t) are: 

m = +-{emr + 0{iemf) , (10) 

A{t) = -^e\Htf + 0{e^{HtY). (H) 

The small parameter in these expansions is e = GA/Stt. We will assume that 
it can be as big as 10~^^ or as small as 10~^^. 

Although deep intuition about the induced stress tensor derives from its 
dependence upon a general metric, we cannot hope to say much about the 
micro-structure of quantum gravity. Nor is this necessary. Screening is a 
phenomenon of cosmological scales, so we need only a model that is accurate 
for spatially flat, homogeneous and isotropic metrics (|^). The isometrics of 
these metrics imply that only Too[6] and Tjj[6] can be non-zero, and that both 
are functions of time alone. We shall parameterize them in the usual way as 
an induced energy density p{t) and an induced pressure p{t): 

foo[6](t) = Pit) , f,,[b]{t) = g,,p{t) . (12) 

We first show that p{t) and p{t) can be chosen so as to support any evolution 
for b{t). Then we show that it is always possible to choose the scalar field 
4>[g] and its potential V^(0) to give the desired energy density and pressure. 
The non-trivial components of the effective field equations are: 

36^ = 3H'^ + 87tGp, (13) 
-26 -36^ = -3H'^ + 8nGp. (14) 

Although one usually regards these as equations for b(t) in terms of p{t) and 
p{t), we can take the converse view: 

Pit) = ^ {3b\t) - , (15) 

Pit) = ^ {-2bit) - 3b\t) + 3H') . (16) 

The physical import of these equations is that one can find pit) and pit) so 
as to support any evolution for bit). 

In a homogeneous and isotropic background, the energy density and pres- 
sure of a scalar (j)[b] are: 

,^\(§)\n$u (IT) 
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P = o -n^)- (18) 



Combining this with (|T5|-0) we see find that an arbitrary evolution b{t) can 
be supported by making the following choices for the scalar and its potential: 




SvrG 
1 



(-26) , (19) 



V = ib + 3b^ - 3H^) . (20) 

Given an explicit function b{t) one constructs V{(f)) by solving the differential 
equation (|l^) for the scalar as an explicit function of time, call it (plbjit) = 
f{t). We then invert this relation to express time as a function of (p, t = 
/~^(0). The potential V{(f)) is found by evaluating relation (^) for t = 

The construction is completed by giving a functional of the metric which 
agrees with 0(t) for the particular choice of b{t). There are many solutions. 
Perhaps the simplest is obtained from P[g]{t,x), the invariant volume of the 
past light cone of the point (t,x). For a spatially fiat, homogeneous and 
isotropic metric, this is a monotonically increasing function of the co-moving 
time t: 

p[b]{t) = J^dt'e^'^'"^ ^'dfe-'^'"^^' , (21) 

so we can invert the relation. Suppose that for the specific function b{t) we 
get P[b](t) = vr(t). Then time is t = 7i~^{P), and the scalar for a general 
metric could be taken as: 

<P[g] = f (n-\P[g])) . (22) 



Since all this can be done for any function b{t), we lose nothing by assuming 
that the induced stress tensor has the scalar form (^. 

In the preceding discussion we inverted the proper order of things to 
show that the scalar stress tensor (^ can describe any homogeneous and 
isotropic geometry. That point having been made, we can return to the 
usual dynamical problem of inferring g^j^^ from T^jj[g]. For a homogeneous 
and isotropic universe the non-trivial equations are: 

3b' = ^H' + 8^gIU^X + V{<P)} , (23) 
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-2b -3b'' = -3H' + 8nG{^(^^^ "^(0)^ (24) 

In dynamical terms, equation (p3|) is actually a constraint. If it is true 
initially then time evolution and energy conservation conspire to keep it true. 



The dynamical equation for b{t) could be taken to be (p4D, but it is more 
convenient to add (|23|) : 

Note that the scalar potential has dropped out. It therefore follows that a 
model is specified by giving the scalar (f)[g] as a functional of the metric. 

If the potential is desired it can be determined by the round-about process 
of first solving (pSj) for b{t) and substituting to find (f)[b]{t). One then inverts 
to express t as a function of (p, and finally substitutes into the constraint 

my- 

2 



v($) 




This turns out to be much easier than it might seem. In Section 7 we will 
obtain analytic expressions for in the perturbative regime and in the 

regime of asymptotically late times. We will also carry out the process nu- 
merically over the full range of evolution. 

It remains to work out during the perturbative regime. Using relation 
d^) we can express the second time derivative of b{t) in terms of the functions 
D{t) and A{t): 




b{t) = H 2D{t) + tD{t) + - ^ - - ^ . (27) 




It is easy to see from the perturbative expansions (|10| ) and ( pAf ) that only the 
third term matters at any time during the perturbative regime . Compar- 
ing with expression (pSf) we infer: 



(28) 
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Making an arbitrary choice of sign and using the fact that only the first term 
in the expansion of A{t) matters, we obtain the following formula for the 
scalar during the perturbative regime: 

^ -^i= ln[l + A(t)] , (29) 



1 



In 



where the small parameter is e = GA/(37r). 



(30) 



3 General principles 

In the previous section we saw that, on cosmological scales, the induced stress 
tensor can be taken as that of a scalar field, (f)[g] which is itself a non-local 
functional of the metric. We saw further that it is really only necessary to give 
this functional, since the associated potential is determined by conservation. 
In fact we actually require only the restriction (j)[b] of this functional to a flat, 
homogeneous and isotropic geometry. However, powerful constraints exist on 
how the scalar can depend upon a general metric. The purpose of this section 
is to state these constraints. 

1. Causality 

We are actually going to guess the induced stress tensor but, were we 
to compute it, T^^[g]{x) would come from Schwinger's effective action for 
expectation values JTO], |ll|], not from the more common, "in"- "out" effective 
action. One consequence is that T^^[g]{x) — and hence also (f)[g]{x) — can 
only depend upon the metric at points which lie within the past light 
cone of x^. It is worth noting that the effective field equations for "in"- 
"out" matrix elements must be symmetric: if they depend upon a field at 
x'^ — Ax^ then they must also depend upon fields at x'^ + Ax^. This is avoided 
in Schwinger's method because his effective action really depends upon two 
fields, one the background during forward time evolution and the other the 
background during evolution back to the initial state. The effective field 
equations are obtained by varying with respect to (either) one of these fields 
and then setting them equal. This is what breaks the forward-backwards 
symmetry of the "in"- "out" field equations. Causality arises because the 
forward and backward evolutions interfere destructively outside the past light 
cone. 
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2. General Coordinate Invariance 

Because the dynamics of quantum general relativity are general coordi- 
nate invariant, non-invariance can only enter the effective field equations from 
the gauge in which the initial state was specified. In other words, (f)[g] must 
be invariant up to surface terms. This means that it must consist of a co- 
variant local part plus non-local operators, such as the retarded propagator, 
acting on local functions of the Riemann tensor and its covariant derivatives. 
Although one can form many covariants from the Riemann tensor and its 
derivatives, only a few are distinct in a homogeneous and isotropic universe. 
To see this, note first that the spacetime is conformally flat. This means the 
Weyl tensor vanishes and one can express the Riemann tensor in terms of 
the Ricci tensor: 

Rp(Tp,v 2 pli^TV Q/JyCrRi/p ~l~ QavRpp, QvpR^a^ g (^9pp-9cn' 9pv9<Tp^ R ■ 

(31) 

Although we will express (j)[g] as an invariant functional of a general metric, 
we need not worry about the distinction between ansatze which agree for a 
spatially flat, homogeneous and isotropic universe. 

3. Stability of the initial value problem, 

Since the quantum field theoretic problem was well posed given only the 
initial state wavefunctional, it must be that the associated effective field 
equations can be evolved forward from t — knowing only the metric and 
its first time derivative. This hmits the local terms in Tnu[9] to those which 
contain at most second time derivatives of the metric. Since the induced 
stress tensor contains derivatives of the effective scalar, the local part of (f)[g] 
can have at most first derivatives. Stability also imposes requirements on 
non-local terms which are differentiated or which can give local terms by 
partial integration. 

4- Non-locality 

A universe which is initially inflating will continue to inflate unless 
is a non-local functional of the metric. To see this, assume the converse. 
Using the previous principle we can then constrain T^z/[fi'] to consist of sec- 
ond rank functions of the Riemann tensor. Note that in a locally de Sitter 
geometry the Riemann tensor can be written in terms of the Ricci scalar. 
One way of expressing this is by saying that the following tensor vanishes: 

^CT/ii/ = Rpafiu Y2 i9pf^9cn/ 9pi>9crtj) R • ("^2) 
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This fact can be used to write any second rank tensor function of the Riemann 
tensor as a term which vanishes in a locally de Sitter geometry plus a function 
of the Ricci scalar times the metric. For example, consider the partially 
contracted product of two Riemann tensors: 

K^'R.ap, = V^^^^V^^p, + \r V\^, + Y^B'g^. . (33) 

But the initial condition of our problem is a locally de Sitter geometry. There- 
fore all terms of the first type vanish, and terms of the second type simply 
renormahze the cosmological constant. We have already defined our cosmo- 
logical constant to absorb any terms of the second type, so local effective 
field equations would have i?^,^ = t^g^^y as a solution for all time. Since 
we can actually see infiation begin to slow in perturbation theory it follows 
that non-local terms must be responsible. The same argument works later 
on, even after the effective Hubble constant has been substantially reduced: 
the effect would stop without non-local terms. So the important part of the 
induced stress tensor must be non-local. 

5. Dimensional Analysis 

The induced stress tensor has the dimensions of length"*^, so (f)[g] goes like 
length"^. This seemingly trivial fact conceals a surprisingly powerful con- 
straint. The most important quantities from which (f)[g] can be constructed 
have the following dimensionalities: 

G ~ length^ , Rp^^^ ~ length"^ , 

A ~ length"^ , ^ ~ length^ , (34) 

where □ is the scalar d'Alembertian: 

n = -^d^ {g^'V^ dy) . (35) 

Note that the dimensionless quantity GA is less than about 10~^^, even for 
GUT scale inflation. Further, the curvature is guaranteed to be of order A 
during the perturbative regime, and it had better be considerably smaller at 
late times. This means that terms with too many powers of G are likely to 
be negligible unless they accompany non-local growth factors such as !/□. 

6. The Flat Space Limit 

When A = we know that the vacuum is stable, so the "in" -"out" matrix 
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element of an operator agrees with its expectation value. We also know that 
"in"- "out" amplitudes with C. loop amplitudes contain at most I infrared 
logarithms |^. This means that the most infrared singular term which can 
remain in the A = "in" -"out" effective action has the form: 

TflatM ~ / d^x^ R [\n{n)fR. (36) 

Terms whose field dependence is more singular must possess positive powers 
of A, for example: 

A j d^x^/^ R^R . (37) 
And we must of course avoid inverse powers of A. 



4 The physics of screening 

The most important constraint on the functional (l)[g\ is that it should cor- 
rectly reflect the physics of screening. Of course choosing the scalar so that 
its perturbative restriction agrees with (^) automatically enforces this dur- 
ing the perturbative regime, so any additional information must come from 
understanding the mechanism of screening for an arbitrary homogeneous and 
isotropic background. That is the purpose of this section. Our procedure is 
to work flrst in the classical background, where results can be checked against 
perturbation theory, and then generalize. We begin by giving a simple deriva- 
tion of the phenomenon of superadiabatic ampliflcation ^, whereby the 
0-point energy of infrared graviton modes is vastly enhanced over the familiar 
of flat space. The next step is to work out the Newtonian approximation 
for the gravitational self-interaction of this 0-point energy. Comparison with 
the known perturbative result indicates how to include relativistic effects. 
Then the analysis is generalized for an arbitrary homogeneous and isotropic 
background. 

Let us recall some facts about the classical background: 

c^^ciass = -dt^ + e^^^dx ■ dx . (38) 

Because it is not physically sensible to assume that coherent inflation com- 
mences over a region of more than about one Hubble volume, we work on the 
manifold x 3?, where each of the coordinate radii is H~^. The 3- volume 
of this manifold is flnite but expands exponentially: 

V{t) = ff-^e^^* . (39) 
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By setting ds^y^^ = we find the world line of a light ray which passes 
through X — 0, directed along the unit vector r at t — to: 

r(t) = 4 (e"''*° - e-^*) . (40) 

Multiplying by exp{Ht) and taking the norm gives the physical distance from 
the origin along the surface of simultaneity: 

e^*||r(i)|| = 4 (l - e-^(*°-*)) . (41) 

From this we infer the existence of a causal horizon of physical distance 
beyond which even a signal traveling at the speed of light can never reach. 
We can also compute the invariant 4- volume of the past hght cone from the 
point (t, x) — {to, 0) to the surface of simultaneity at t — 0: 

;*« 4 

class 1 ^-0 



(to) = f ' dt e'""' X -TTimr , (42) 
Jo 6 

" — ' (wto - ^ + ie-"" - ^e-™« + ie-«"«) . (43) 



Finally, note that the coordinate transformation 77 = —H~^ exp^—Ht) makes 
the classical background conformal to flat space: 

dsl^,, = {-dv' + dx ■ dx) , (44) 

where the conformal factor is = —l/{Hri). Note that the surface of si- 
multaneity at t = corresponds to 7] = —1/H, and that the infinite future 
corresponds to — > 0~. 

The next step is the kinematics of free gravitons. Graviton modes are 
described by a polarization and by a co-moving wave number of the form 

— * " 

k = 27rHn, where n is a 3-tuple of integers. The integral approximation to a 
mode sum is: 

2 X: / {27rHn) ^ 2 J d'n f {27rHn) = 2 J f{k) , (45) 

where the infrared cutoff is at k = \\k\\ = H. Since physical distances expand 
by exp{Ht), physical wave numbers redshift by exp{—Ht). We are most 



12 



interested in infrared modes, defined as those whose physical wave lengths 
have expanded beyond the causal horizon: 

Infrared ^ H <k< He"' . (46) 

We shall refer to the higher modes as "ultraviolet." 

Now consider the dynamics of free gravitons. For any homogeneous and 
isotropic geometry, these are the same as those of a massless, minimally 
coupled scalar Suppose we call such a field ipijj^x). In the classical 
background its Lagrange density is: 

C = hl''(ij'' -Vij-Vij) , (47) 

where a prime denotes differentiation with respect to the conformal time rj. 
The mode coordinates are obtained by taking the spatial Fourier transform 
and multiplying by a factor of H: 

q^(7]) = H j e^^-^V'(r/, x) . (48) 

These variables allow us to recognize the Lagrangian as a sum of independent 
harmonic oscillators: 

L^jd'^ = \h^' E {^U - k% g,^) . (49) 

k 

Since g_g = gi we can treat this system as if there were a single real mode 

for each wave number k. 

It is straightforward to express the mode coordinate and its conjugate 
momentum in terms of creation and annihilation operators. The Heisenberg 
equation of motion is: 

0. (50) 
It follows that the negative frequency mode solution is: 

u{r^,k) = ^{l-^\e~''\ (51) 



'2k \ krj 
and we define its time derivative as: 



M'(r/, k) = -ikn~\{ri, k) = -ik^=e-'^~^ . (52) 
13 



The Wronskian formed from u{ri,k) and v{r],k) is constant in consequence 
of the equation of motion, and with our normahzation its value is: 

u{ri, k)v*{ri, k) + m*(?7, k)v{ri, k) = — . (53) 
We can express g'^(?7) as a hnear combination of u{rj, k) and u*{ri, k): 

q^ivi) =u{ri,k)a^ + u*{r],k)a\. (54) 
The conjugate momentum is: 

p^^i^l) = Hn^'^ir^), (55) 
= —ikHv{r], k)ai: + ikHv{r], k)al . (56) 

Requiring that 5'^(?7) commute canonically with p^(?7) and making use of the 
Wronskian determines the commutator of ar and at to be: 

The "Hamiltonian" which generates the conformal time evolution of mode 
k is: 

*E = Jp4^*^"\^ (58) 

This is just a harmonic oscillator with frequency k and mass HVP'. Because 
the mass is time dependent there are no stationary states but one can of 
course compute the expectation value of in the presence of some state. 
In the far ultraviolet curvature is obviously a small effect, so we may assume 
the flat space vacuum: 

agio) = . (59) 

Note that Heisenberg states do not evolve, and that the operator ag was 
constructed to be time independent. Hence condition (|^ ) persists, even 
after the originally ultraviolet mode has red shifted to the infrared. 

The expectation value of is simple to take in the presence of this state: 

mi^m = ^xH'kMv,k)v*{r^,k)x^ 

+]-Hk'^n^ X u{ri, k)u*{ri, k) x ^ , (60) 
2 H 
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(61) 
(62) 

Now exploit the relation between co-moving time and conformal time to 
relate the co-moving Hamiltonian to the conformal one: 

if* = = in-i^ = (63) 
dt dri ^ ^ 

— * 

It follows that the physical energy in mode k at co-moving time t is: 

The first term is just the familiar 0-point energy, appropriately red shifted. 
One way to understand the second term is that virtual gravitons whose phys- 
ical wave lengths exceed the Hubble radius cannot recombine; they are pulled 
apart by the expansion of spacetime. The energy of any one such graviton 
redshifts, but there are so many produced that the total energy contributed 
by each infrared mode actually increases. 

Since there are an infinite number of ultraviolet modes, the total 0-point 
energy diverges. This is not consistent with the assumption that the back- 
ground is initially undergoing inflation with Hubble constant H. To make 
the assumption consistent we must subtract the original 0-point energy by 
normal ordering. However, this has only a minuscule effect on the superadi- 
abatically amplified 0-point energy of the infrared modes. We can use simple 
Newtonian ideas to obtain a crude estimate of the energy density induced by 
their self-gravitation. 

To obtain the energy density of mode k we divide by the 3- volume: 

P.-(^.--^^e--*)^nt) = ^e--*. (65) 

Although this red shifts to zero, it does so more slowly than pure radia- 
tion. One consequence is that the associated Newtonian potential remains 
constant: 

- e-^^*/cV,— 47rGp,^ ^ ^ . (66) 



15 



The total Newtonian potential from all infrared modes is accordingly: 

Vi^ = ^im dkk^Vu = Ht . (67) 

This combines with the total 0-point energy density of infrared modes: 

1 rHcMHt] ,,,2/74 
^^^ = ^7^ ^^^P^ = ^' (68) 

to produce an increasing gravitational interaction energy density: 

PNcwt(t) = fm ■ Pm = — ^^Ht . (69) 

This is down from piR by a factor of the small number GA/(37r) < 10^^^, 
but its time dependence eventually makes it the more important effect. 

The Newtonian estimate we have just obtained compares fairly well with 
the exact result of the lowest non-trivial order in perturbation theory [|I[]: 

Pit) = {{Hty + Om} + 0{G') . (70) 

The extra factor of Ht derives from the inclusion of four relativistic effects 
which were omitted in the Newtonian approximation: 

1. There is a 0-point pressure in addition to the 0-point energy density. 

2. The Newtonian potential is not the only gravitational field. 

3. The various gravitational potentials can carry momentum. 

4. The gravitational interaction is not linear. 

Although one must really do the quantum field theoretic calculation to get 
the right answer, it is simple enough to indicate how each of these effects 
modifies the Newtonian estimate. 

If we assume that the 0-point stress-energy of each mode is separately 
conserved then the 0-point pressure works out to be: 
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Like the 0-point energy density, the total 0-point pressure is only a constant 
— and very small — fraction of the pressure in the cosmological constant. 
However, unlike the 0-point energy, the 0-point pressure serves as the source 
for a gravitational potential whose homogeneous equation of motion is that 
of a massless, minimally coupled scalar |jl2|. It is straightforward to compute 



the retarded Green's function for this potential: 

G iv, x; v', ^) = -2d{v -v')J-0^, e'<'~''hm [u%v, k) u{^\ k)] , (72) 
= JAM. I ^'2^' ^ ( Ary - Ax) + H^e {At] - Ax) | , (73) 

where Arj = t] — rj' and Ax = ||x — x'\\. The first term is just like its fiat 
space cognate: the only contribution comes from sources on the actual light 
cone, so there is no growth for a constant source. However, the second term 
superposes over the entire past light cone, whose invariant volume (^3]) grows 
like Ht in the classical background. 

The third relativistic effect means that we should not view the precip- 
itating event as the creation of two gravitons with opposite 3-momenta, 
whose stress energy then induces a gravitational potential containing zero 
3-momentum. What really happens is the creation of a graviton with 3- 
momentum ki and another with 3-momentum k2, which together induce a 
potential with 3-momentum — (fci + ^2)- Although this still leaves two mode 
sums, they are not cleanly split between a 0-point stress energy and a poten- 
tial term, as was the case for our Newtonian estimate. 

The final relativistic effect means that we must include a bewildering 
variety of interactions where the potential scatters off one of the gravitons, 
or where it interacts with itself. This is one of the things that makes the 
quantum field theoretic calculation so difficult. We can nonetheless say that 
the effect is still due to the gravitational interaction between virtual gravitons 
which are pulled apart by inflation. Superadiabatic ampliflcation can still be 
used to estimate the rate at which these gravitons are created and the stress 
energy which they carry. And it is generally the case that one factor of Ht 
derives from one of the two mode sums while the other factor of Ht comes 
from an integration over interaction times. 

Correcting the Newtonian estimate for the induced energy density is not 
quite the end. Just as a pressure is associated with p^, so there is a pressure 
associated with p{t). We can flnd it by using energy conservation, which 
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reads as follows in the classical background (^: 



p = -m{p + p). (74) 

Combined with ( |7DD this implies that most rapidly growing part of the in- 
duced pressure is exactly opposite that of the energy density: 

pit) = +^ {{Htf + Om] + O(G^) . (75) 

It follows that the interaction between infrared gravitons acts to screen in- 
flation by an amount that becomes non-perturbatively large at late times. 
At this point it is useful to recall some standard facts about inflation 



I3| , |T^ in order to form a proper impression both of the effect's magnitude 
and of the time scale over which it acts. What is usually termed, the "scale 
of inflation," is the mass M defined so that equals the energy density of 
the cosmological constant, A/(87rG'). Since the Planck mass is Mp = G^^/^ 
we have: ^ 

GK = ^'k(^\ . (76) 



\Mp/ 

It is traditional to assume that M is GUT scale, which makes GA ~ 10^^^. 
One sometimes encounters models with scales as low as that of electroweak 
symmetry breaking. (Past that point there is not enough CP violation to 
explain the observed baryon asymmetry.) Inflation on the electroweak scale 
would give GA ~ 10"^''. These numbers mean that the gravitational inter- 
action energy density is a very small faction of M'^ unless Ht is enormous. 

Although there are higher order corrections to the induced energy density 
( [TDI ) and pressure (|75|), it turns out that the lowest order effect becomes 
non-perturbatively strong when these higher terms are still insignificant 0. 
The way this works is that the induced stress tensor serves as a source for 
corrections to the classical background (H). The pressure again engenders an 



extra factor of Ht from the invariant volume of the past light cone [0, and 
this causes the lowest order effect to throttle inflation before higher orders 
can become significant. When the background is expressed in co-moving 
coordinates (0), the function h{t) has the form: 

h{t) = Ht[l + Dit)] + i ln[l + Ait)] , (77) 
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where we recall from ( PI1| - |TT| ) the perturbative expansions of D{t) and A{t) 
for small e = G'A/(37r) and large Ht [^]: 

m = +-{eHtf + 0{{eHtf) , (78) 

A{t) = -^e\Htf + O {e^Htr) , (79) 

Perturbation theory breaks down when A{t) approaches —1, at which time 
only the first term in the expansion of A{t) is relevant and no term in D{t) 
is significant p. 

The function A{t) is of great importance because we saw, at the close of 
Section 2, that it gives the scalar during the perturbative regime: 

0(t) = --i=ln[l + A(t)] . (80) 
V ottG 

Summarizing and abstracting the preceding analysis, we may conclude that 
A{t) derives, during the perturbative regime, from acting the retarded scalar 
Green's function, □ , on a source that grows like {HtY: 

A = -SttG f'dt'e'^"'' f dt"e^"'"SouTce{t") = StigI- (Source) , (81) 
Jo Jo □ 

Source = 172 ^{Htf + ... . (82) 

This source is the stress energy induced by the gravitational interaction be- 
tween gravitons produced by superadiabatic amplification. It consists of 
a variety of stress energy-potential and potential-potential terms, each of 
which contains two mode sums over the 3-momenta of the gravitons. One 
factor of Ht is attributable to one of these mode sums, the other comes from 
the superposition over interaction times in a potential. Part of the task of 
generalization is straightforward since the potentials are obtained from the 
retarded Green's function, which can be defined for any homogeneous and 
isotropic background: 

!/ = _/* rft'e-3^(*') /*' dt"e''^'"^f{t") . (83) 

n JO Jo 

What remains is to understand superadiabatic amplification on a general 
homogeneous and isotropic background. 
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We begin by comparing the co-moving element with the conformal one: 

ds'^ = -dt^ + e^^^^'^dx ■dx = Vl? (-dr{^ + dx ■ dx) , (84) 

to infer general relations for the conformal factor and the coordinate trans- 
formation: 

n{r]) = exp[b{t)] dt = Vtdr] . (85) 

Let us denote hj rj = rji the image of t = 0. Although we cannot be explicit 
for a general b{t), any superluminally expanding spacetime will have rji < 0. 
We can also assume that the approach to the infinite future is ?7 — 0~, as 
before. The manifold is still x 3?, with each of the coordinate radii equal 
to H^^. We therefore conclude that the 3-volume and the graviton mode 
sum are: 

V{t) = if-^e^^W , (86) 
2Y.f{27rH\\n\\) ^ J-^ J dkk'fik) . (87) 

n 

It turns out that there is always a factor of V^^ associated with each mode 
sum so that the factors of cancel and any physical dependence upon the 
range of co-moving coordinates must come from the limits of integration. 

Dynamical graviton modes still obey the equation of motion of a massless, 
minimally coupled scalar ||^ . The mode coordinates for such a scalar are still 
obtained by spatial Fourier transforming according to (^8]), and the formula 
for the Lagrangian is unchanged from (^91) , provided the general conformal 
factor is understood. What changes is the Heisenberg equations of motion: 

+ 2§?;^ + k\ = . (88) 

Redefining the field variable as Q{ri) = VLqj:{vi) gives the following suggestive 
equation: 

Q"+(fc'-^)Q = o. (89) 

There are two regimes in which good approximate solutions can be found: 
the far ultraviolet, where k"^ ^ and the far infrared where k"^ ^ Vt" /VL. 

Before going on to study these regimes we should note the important fact 
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that they can be given an invariant specification. In view of the relation (53) 
between conformal and co-moving coordinates we can write: 

^" ^e'^e' = e'' (b + 2b') = -e^'R . (90) 



Q dt dt ^ ^6 

The two regimes are therefore characterized by how the physical (i.e., red 
shifted) momentum compares with the square root of the Ricci scalar. In 
view of our definition ( ^6l) for the classical background, it is reasonable to 
make the following general definition for "infrared" gravitons: 

Infrared ^ e-^'^'^H^ < e-^'^'^k' < ^R{t) . (91) 

As before, we refer to the higher modes as "ultraviolet." 

In the ultraviolet regime the Ricci scalar term is a perturbation: 

Q" + k'Q = . (92) 
We can obtain the mode functions by iterating those of fiat space: 

Q^^ir^, k) = ^e-^^^ + dr] - sin [k{r^ - f])] QUv. k) . (93) 



This obviously results in a series in powers of l/fc, the successive terms of 
which are less and less significant for large, negative rj. 
The mode functions associated with Quviv^ k) are: 

u{7],k) = n'\7])Q^,{r],k) , (94) 

viv.k) = '-[niv)Q'^^i^,k)~n'iv)QnAv,k)] . (95) 

The associated Wronskian: 

m(?7, k)v*{ri, k) + u*{rj, k)v{rj, k) = 

[QUv, k)QUv, k) - QUV, k)Qt{v, k)] . (96) 



is constant in consequence of the equation of motion (|89D. Since Qmiv^k) 
approaches the fiat space mode functions at early times, we see that the 
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constant is the same as before. The subsequent operator expansions 
and commutation relations: 



u{r], k)a^ + u*{7], k)ai , 
—ikHv{r), k)aj^ + ikHv{r], A;)at , 



1 

H 



(97) 
(98) 

(99) 



are the same as for the classical background. And the expectation value of 
the Hamiltonian in the presence of the state annihilated by is: 



1 



Quv Q^^^ 



k Qy^yQuV ■ 



(100) 



Although this formula is correct for all times, obtaining the leading behavior 
for small r) requires that we develop an infrared expansion for the mode 
functions. 

In the infrared regime the momentum term is a perturbation: 



At zeroth order the two independent solutions are: 



dr] 



;ioi) 



(102) 



For superluminal expansion the limit 77 — 0~ carries Qioii]) to infinity, while 
Q2o{v) goes to zero as ~ rjQ^^. The limit f] — > —00 has the opposite effect: 
Qioiv) approaches zero and Q2o{v) go^s to infinity as ~ rjQ^^. 

In view of the limiting forms for QiqIt]) and Q2o{v) it is the advanced 
Green's function: 



Gadv{r],v) = ^iv - V) [Qioiv)Q2oiv) - Q2oiv)Qioiv)] , 



(103) 



which gives a reasonable equation to iterate for the full solution based on 

Q20- 

Q2{v,k) = Q2o{v) - f drjG^,{rj,rj)Q2{fJ,k) . (104) 

J—oo 
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The resulting series is: 



Q2 = E(-^')'G'ldv ■ Q20 , (105) 
where the £-th "power" of the Green's function denotes the £-fold integration: 
G'adv-Q2o=/ c??7i Gadvl??,??!) ■ ■ ■ / drji G i,d^,{r]i_ i, 7] i)Q 20 {ve) (106) 

For asymptotically small rj this £-th power goes to zero like ~ 7f^VL~^. The 
expansion ( |105| ) is therefore in terms which are less and less significant at 
late times. 

The square of Qwiji) is not integrable at = 0. To obtain the full solution 
based on Qio we must therefore begin iterating with a Green's function which 
is intermediate between advanced and retarded: 

GUv,v)=d{v-v)Q2o{v)QiM + G{v-ri)QMQ2M ■ (107) 

The result is a series whose terms have the asymptotic behavior ~ rf^VL. For 
high enough £ the integral of such a term times Qiq no longer converges at 
large rj. At this point we must continue the iteration using the advanced 
Green's function. The full series is: 

N oo 

Qi = T^i-k^GL ■Qio+ E i-k'YGi^.'' ■ GZ, ■ Qio , (108) 

e=o e=N+i 

where N is the order of perturbation theory at which we must change Green's 
functions. As with the expansion for Q2, the small rj behavior of each suc- 
cessive term is weaker than that of its predecessor by a factor of rj'^. 

Since Qi(?7, k) and Q2(^, k) are independent solutions of the same linear, 
second order differential equation as QuviVy k), we can find complex numbers 
a and (3 to enforce the condition: 

QUv.k) = aQi{r],k)+PQ2{v,k) , (109) 
= QuAv,k). (110) 

Although Qij-{ri,k) = QuviVjk) over the full range of 77, we can estimate a 
by matching the zeroth order solutions at the boundary between ultraviolet 
and infrared. This gives the relation: 

a^^e-''""', (111) 
V2k 
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where t]i is the conformal time at which 2k'^ = Q'l/Qi. 

Our estimate ( p. 1 1|) can be checked whenever exphcit solutions exist. For 
example, suppose the scale factor is a power law which obeys the same initial 
conditions (&(0) = and 6(0) = H) as the classical background: 

em=(^l + Bj . (112) 

(We assume p > 1 to make the expansion superluminal.) The conformal 
factor is: 

p 

n=(^^y~\ (113) 

where rjQ = —p/{p — 1) and the properly normalized solution is: 

Q(r^, k) = ^v^e-'f ('^+^)if(^)(fcr^) , u = ^- (j^) ■ im 

The infrared solutions Qi{ri,k) and Q2{ri,k) give the two Bessel functions 
from which the Hankel function is constructed: 

aQi{'n,k) = ^csc(i/7r)e-^('^+^)^0F?yj_^(fcr/) , (115) 

PQ2{v,k) = -^csc(z/7r)e5('^+5)-0F^J,(/£r^) . (116) 

By making the first terms of the respective series expansions agree with 
Qwiv) Q^oiv) W6 infer the following expressions for a and /5: 

r(^ /z^\ ^ , (nr) 

V2^\kVo) V32 8 ; 72^' ^ ^ 



'2% csc(z/7r) ( ikriQ\^ ^ 

Tr(z/ + i) \ 

csc(z^/i 1 ^y'^-i nr' 

'Iviy/Ti— V . (lloj 



rz/ + l V32 8 7 J2k 



Except for powers very close to the superluminal bound (p = 1) the index 
V is of order one and we see that a agrees with estimate (|111|). Unless u 
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happens to be very close to an integer or half integer we also see that a and 
/3 are comparable. 

Our series expansions imply the following relations for late times: 

QUv,k) ^ an{r]), (119) 

Q'(n) 

Q[r{v,k)--f^QUv,k) ^ av^iv). (120) 
i2(r/) 



Substitution in ( |100|) gives the late time behavior of the (conformal) energy 
in mode k: 

miivm - laa*eQ\v) ~ , (121) 

where -Ri is the value of the Ricci scalar at the time mode k crosses from the 
ultraviolet regime to the infrared. The form of Ri{k) depends upon b{t). For 
general power law inflation ( |112| ) the Ricci scalar is: 

m ^ 6 (m + 2i^t)) = [j^)' [-1 + 12) , (122) 
and we find: 

2p-l\'- (H\— 



Taking p to infinity recovers the de Sitter result, Ri{k) = 12H^. We therefore 
expect that Ri{k) is generally a slowly decreasing or constant function of k. 

It is now straightforward to generalize what was done in classical back- 
ground. The co-moving energy is down by a factor of Q so it approaches: 

^ ^ laa*ee'^'^ ~ ^e^W . (124) 
^2 48k ^ ^ 

Dividing by the 3-volume gives the energy density: 

p, - iaa*A;^if^e-^« ~ ^e-^« . (125) 

— * 

If the stress energy in mode k is separately conserved it is straightforward to 
infer the pressure: 

Pk = -^KPk+Pk) =^ Pk = -\pk- (126) 
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We can even parallel our previous estimate for the Newtonian potential of 
mode k: 

-e-^WfcV.— 4vrGp,^ ^ ^ = '^TrGaa* ^ -^^^^ . (127) 

Of course all these results reduce to those obtained for the classical back- 
ground when we set Ri = 12H^. The difference is that we can now recognize 
the Ricci scalar as the source of the stress energy and the gravitational po- 
tentials whose various interactions lead to screening. 

It remains to see what becomes of the two mode sums. Although we have 
already explained that the Newtonian estimate misses some essential features 
of the fully relativistic result, it can still be used to get a rough idea. In the 
Newtonian approximation the total infrared energy density and potential are 
the following mode sums: 

PiR = ^773 / dkepk-T^^'^'^^ dkkRiik), (128) 
TT^i/'^ Jh 487r^ Jh 

1 rK(t) n i-K(t) 

- (129) 

where K'^{t) = i?(t) exp[26(t)]/12. The point which generalizes about this 
is that, for the dominant terms, one of the mode sums goes like exp[2b]kdk, 
whereas the other has the form dk/k. 
The relation: 

k' = ^R,e''^ , (130) 

can be used to convert these mode sums into integrations over time. For 
example, the general power law inflation (|112|) gives: 

kdk = ( I —Rie^Hidti . (131) 

\ p ) 12 

Except for very slow inflation (p ~ 1) the time dependence of the Ricci scalar 
is negligible compared with that of the scale factor. We can therefore write: 

kdk ^ —Rie^Hidti . (132) 
12 ^ ' 

In the same approximation the total infrared energy density becomes: 

p,^ ~ ^-e-2^W /* dt'h{t')e''^''m\t') ^ . (133) 

^ " 5767r2 Jo ^ ' ^ ' 11527r2 ^ ^ 
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Hence we conclude that the exp[2b]kdk mode sum engenders an extra factor 
of the Ricci scalar. 

The dk/k mode sum is more subtle. Of course we can convert to an 
integration over time as before: 

<PiR^-^ f dt'h{t')R{t') . (134) 

i/TT JO 

However, further progress requires that we neglect h relative to 1? . For the 
general power law ( |112| ) the ratio h/lP' is — 1/p, which suggests that the 
approximation is also valid for rapid inflation. We can use this to introduce 
a second integration, and write the result in terms of □ : 

m ~ -- fdt'h^t') , (135) 



TT JO 

3G 



— / rft'e-3''(*') / dt"e^'^''%\t") (136) 

TT JO JO 

^ ' [R^) . (137) 



487r □ 



.-1 



Our conclusion is accordingly that the dk/k mode sum goes to □ R. 

It would be a mistake to take the Newtonian approximation too seriously. 
For example, although the source term of equation (|HT| ) contains two mode 
sums, they are not associated one with a 0-point stress energy and the other 
with a potential. Nor must the source have the Newtonian form of a stress 
energy times a potential; one can also have products of derivatives of poten- 
tials, and there are terms in which the two virtual gravitons interact with 
the potential at different times. What can reasonably be concluded is that 
the source term of equation (|8TD involves two factors of D"^ acting in some 
order on five Ricci scalars. 



5 Ansatze for the scalar 

The purpose of this section is to discuss possibilities for the scalar (j)[g] that 
are consistent with the results of the last three sections. We begin by re- 
viewing these results, then we tabulate the 33 candidate invariants which 
can be constructed from the retarded scalar Green's function and the 
Ricci scalar R. The section closes with a discussion of the possibilities for 
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involving other Green's functions, and for replacing some of the five Ricci 
scalars by perturbatively equivalent factors of 4A. 



From Section 3 we learned that F[g\{x) must be a dimensionless invariant 
which depends causally upon the fields within the past lightcone of x^. The 
requirement of a stable initial value problem implies that any factors of the 
Riemann tensor must be protected by at least two inverse derivatives, with 
derivatives of the curvature protected by correspondingly more inverse deriva- 
tives. This is consistent with the fact that the important part of F[g\{x) must 
be non-local. However, correspondence with the known results for A = puts 
severe restrictions on how many inverse derivatives can appear. To be pre- 
cise, the induced stress tensor must either vanish with A or else it must 
not be more infrared singular than 9^ ln(9^) when subjected to the derivative 
expansion for weak fields in a fiat space background. 

Section 4 used the physics of our mechanism to show that F[g\ has the 
form of the scalar retarded Green's function D"^ acting on a source con- 
structed from two more factors of C"^ and five factors of the Ricci scalar 
R. This reduces the problem to combinatorics. There are 21 ways of placing 
the factors of R when the two retarded Green's functions act in series, and 
there are 12 distinct placements when they act in parallel. Table 1 lists the 
33 possibilities. If one requires that at least one factor of R must reside at 
each of the three locations then only ten terms remain. 

All of the terms in Table 1 are manifestly invariant, causal and non-local. 
Note also that various the factors of R are protected by enough factors of 
so as not to jeopardize the stability of the initial value problem. There is 
no problem with the fiat space limit since the source terms behave generically 
like in the weak field derivative expansion. This means that the scalar — 

^In fact F{t) is just minus the perturbative coefficient function A{t) of expression (pi]). 





(Ht)^ + sub dominant , 



(139) 
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# 


Candidate Source 


# 


Candidate Source 


la 




lb 


i?5(n-ii) (n-'i) 


2a 


2R^{n-^R{n-n)) 


2b 




3a 


2R^{n'^i{n-^R)) 


3b 


R^{n-'R^) (n-H) 


4a 


2i?3 (n-'R^n-H)) 


4b 


R^iD'^R) (n-^R) 


5a 


2 R^n-^R{n-^R)) 


5b 


R\n-'R^) (n-^i) 


6a 


2i?3 {n-nin-'R')) 


6b 


R\n-'R^) (n-'R) 


7a 


2i?2 (n-ii?3(n-ii)) 


7b 


R{n-'R^) (n-'i) 


8a 


2R^D-'R^{n-'R)) 


8b 


Rin-^R"") in~'R) 


9a 


2R^{n-'R{n-'R^)) 


9b 


R{n-'R^) (n-'R^) 


10a 


2R^ {n-'i{n-'R^)) 


10b 


(c-H) 


11a 


2R{n-^R\n-^i)) 


lib 


(n-^R^) (n-^R) 


12a 


2R{n-^R^{n-^R)) 


12b 




13a 


2i?(C-ii?2(C-ii?2)) 






11a 


2 R (n-\n(n-\R'^)) 






15a 


2R{n-H{n-^R'')) 






16a 


2 {n~^R^n-H)) 






17a 


2 {n'^R\n-'^R)) 






18a 


2 (□-ii?3(n-ii?2)) 






19a 


2 (□-ii?2(n-ii?3)) 






20a 


2 (□-ii?(n-ii?4)) 






21a 


2 (□-H(c-ii?5)) 







Table 1: Candidate sources for F[g]. To obtain the functional F[g] act 
— tI ( ) □"''^ on the source. 
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and hence the induced stress tensor — goes hke d^, which is perfectly allowed. 
Though it would not have been legitimate to impose this as a requirement, it 
is reassuring to note that all of the terms approach a constant for subluminal 
expansion. 

It is important to note that the actual source need not consist of a single 
term from Table 1; it could equally well be a linear combination. Though we 
will not attempt it here, an explicit derivation of the source seems possible in 
two different ways. First, we might re-compute the leading two loop diagrams 
for a general homogeneous and isotropic background, using the methods of 
Section 4 to isolate the dominant terms and to express them in invariant 
form. This would be difficult, but no more so than the original computation. 
The second technique exploits the fact that the terms of Table 1 survive in the 
limit A — »• 0. They must therefore appear in the usual effective act ion, which 
can be invoked to further constrain the allowed combinations. One might 
even be able to pick off the coefficients by subjecting the two loop effective 
action to an expansion in powers of curvature along the lines developed by 
Barvinsky and Vilkovisky ||T5|. 



It remains to discuss two issues, the first of which is why the non-locality 
is confined to inverses of the scalar d'Alembertian: 




dt\ dtj dt\ dt) ^ ' 

Of course we derived its presence from the physics of the mechanism, but 
why do we not need to allow for the possibility of inverses of other differential 
operators? It has been pointed out that the kinematics of free gravitons on a 
homogeneous and isotropic background is governed by □ but there are also 
non-dynamical, constrained modes which possess a different kinetic operator 

Db ^ -e-'^ ( e-'^eA . (141) 



dt 

We will additionally consider the kinetic operator of a massless, conformally 



^This sounds paradoxical — in view of the fact that the "in" vacuum does not evolve 
into the "out" vacuum — but it isn't really. The distinction between the two vacua 
simply means that the imaginary part of the effective action is non-zero (in fact, infrared 
divergent) when evaluated at the classical background. 
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coupled scalar: 



S„„,,a_li?._e-|(e'|e'), (142) 



although there is no dynamical reason to suspect its presence. 

The reason we ignore D^^ and -D^^f is that they act to produce constants 
during the perturbative regime, and they are even less relevant for slower 
expansion. To see this, consider how each of the three inverse differential 
operators acts on a function of time f{t) for an arbitrary homogeneous and 
isotropic background: 

If = - tdt'e-^'^''U'' dt"e'''^'"^f{t"), (143) 
□ Jo Jo 

= -e-^m ['dt'e'^^'U'' dt"e'^'"^f{t"), (144) 
Db Jo Jo 

J—f = -e-'-W /*di'e-^(*') rdi"e2^(*")/(t"). (145) 

i^conf Jo 

In the perturbative regime b{t) — Ht and R{t) — 12H'^. Therefore acting 
each of the three operators on the Ricci scalar gives: 

= -AHt + -- -e"^^* , (146) 
□ 3 3' ^ ^ 
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l-R = -6(l-2e-^* + e-2^*) , (147) 
J—R = -6(l-2e-^* + e-2^*) . (148) 

-^conf 

Except for the beginning of inflation, the operators Db and i^conf are indis- 
tinguishable from constants during the perturbative regime! On the other 
hand, □ grows hke Ht. When we recall that this factor must be at least 
10^ before anything interesting happens, the other operators are negligible 
in comparison. 

The other operators perform even worse during the slower expansion that 
should follow the end of inflation. Suppose we call the transition time tz- It 
is useful to break the double integrations of the operators into periods which 
are before and after t^: 

ft ft' ft, ft' ft ft, ft ft' 

/ dt' / dt" = dt' dt" + / dt' dt" + / dt' dt" . (149) 
Jo Jo Jo Jo Jtz Jo Jtz Jtz 
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Since there are so many more e-foldings before tz than afterwards the largest 
contribution for each of the three operators comes from the first of the three 
integrals to the right of ( |149| ) . Since this integral is restricted to the per- 
turbative regime we can estimate the post inflationary behavior of the three 
operators: 

^ -iHtz , (150) 

ii? ^ _Qg-2f,W+2Mt.) (151) 

Db 

J—R ^ _6e-^W+''(*^) . (152) 

Dconf 

Whereas □ ^ gives an enormous constant, the other operators are only of 
order unity at t = t^, and they fall off thereafter. It follows that these 
operators are irrelevant. 

So much for other operators; the final issue is whether some of the five 
factors of R should be replaced by 4A. Although these terms are distinct 
for a general homogeneous and isotropic background, it is not possible to 
distinguish them during the perturbative regime: 

R — >R = Q(b + 2b'^) — >R = 12H'^. (153) 

For the factors of R which occur furthest back in time it probably makes 
little difference whether or not they are replaced by 4A. The integrals are 
dominated by the inflationary period, during which the two terms agree. On 
the other hand, the factors of R which appear latest might play an important 
role in the post inflationary regime where they are insignificant compared 
with 4A. 

Of course one can consider the physics of screening on a general back- 
ground, as we did in Section 4. This favors R over 4A. The preference is 
very strong for the two factors of R from the 0-point energy of the virtual 
gravitons ( |125D and for the factor from the kdk mode sum ( |132| ). It is some- 
what weaker for the factor of R from the dk/k mode sum ( |137|) , and it is 
very thin for the factor that was introduced to compensate the dimensions 
of the other in the source. 

We can therefore find some justification to consider sources involving four 
or only three factors of the Ricci scalar. There are 15 ways of placing four 
i?'s when the two factors of C"^ act sequentially, and there are 9 distinct 
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# 


Candidate Source 


# 


Candidate Source 


Ic 


/-X 7-vQ 1 1 1 1 1 -1 \ \ 

2 i?'^ (n ^i(n ^1)) 


Id 


R'^{n ^1) (□ ^1) 


2c 


r» 7^9 1—1 7->/ 1 — 1 -1 \ \ 

2 i?^ (□ '■R{n ^1)) 


2d 


^-R) (□ ^1) 


3c 


2 i?^ (□ ^i?)) 


3d 


7-> / 1 7->9\ 1 1 1 \ 

R[n '^R^j (n ^1) 


4c 


r» / 1—1—1 n2/i-n — 

2 i? (□ ^it^(C ^1)) 


4d 


i?(C -R) (□ it!) 


5c 


2i?(n-ii?(n-ii?)) 


5d 


(u-^R^) (n-H) 


6c 


2i?(c-ii(n-ii?2)) 


6d 




7c 


2 (n-ii?3(n-ii)) 






8c 








9c 


2 {n~'R{n~'R^)) 






10c 


2 (□-H(n-ii?3)) 







Table 2: Candidate sources for ^[^f] with only three Ricci scalars. To obtain 
the functional F[g] act — □"^on the source. 



placements for the parallel case. We shall not bother listing them all. Ta- 
ble 2 gives the 10 series and 6 parallel candidates that contain only three 
factors of R. It is worth remarking that the presence of explicit factors of 
A would preclude deriving the induced stress tensor from the usual effec- 
tive action, however, one could still derive the actual source by re-computing 
the dominant two loop contribution in a general homogeneous and isotropic 
background. 



6 Numerical evolution 

Despite the enormous restriction we have obtained on the possible forms for 
the scalar (f)[g], it is still not quite unique. The purpose of this section is to 
develop and to implement a scheme for numerically evolving whatever model 
is defined when one finally settles upon a choice for (f)[g]. Our procedure will 
be to select one of the ansatze of the previous section and then evolve it 
far enough past the end of inflation to infer the asymptotic behavior at late 
times. We shall also obtain explicit results for the number of e-foldings which 
occur before the end of inflation, the rapidity of the transition, and for the 
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equation of state throughout the evolution. It ought to be straightforward 
to generahze the numerical scheme to any of the viable ansatze discussed in 
Section 5. 

Recall from Section 2 that the evolution equation is: 

subject to the initial condition b{0) = and 6(0) = H. This gives b{t). From 
the field effective field equation ( p!^ - |T^) we see that the total energy density 
and pressure, including the contribution from the cosmological constant, are: 

P» = 8^3i^ (155) 

Pto. = (-26 - 3b^) . (186) 

This means that the instantaneous equation of state can be reconstructed as 
follows: 

P^t2^_26(t)_^_ (157) 

Recall as well that the scalar can be written in terms of a dimensionless 
invariant F[g]: 

0b] = -^i^ln(l-FM) . (158) 
V ottG 

Substitution gives the following general evolution equation: 

2\l- F\b]) 

In Section 5 we argued that F[g] has the form of a retarded scalar Green's 
function acting on a source composed of two more retarded scalar Green's 
functions which act in some order on three to five factors of the Ricci scalar. 

Before specializing to a particular ansatz, we should comment on the 
qualitative behavior of all models. During inflation the time derivative of the 
dynamical variable b{t) is a large positive constant. Inflation is brought to 
an end when the inherently negative second derivative becomes large enough 
to force b{t) down to nearly zero. It is clear from (|1 591) that this must occur 
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when F[b]{t) approaches unity. All viable models show the following growth 
for F{t) during the perturbative regime: 

172 

F{t) = —e^{Htf + subdominant , (160) 
o 

where we define the dimensionless parameter e = GA/ (Svr). Our perturbative 
estimate for the number of infiationary e-foldings is accordingly 

'^--(m) (ga) "(tIoos) (if) • 'i"' 

where Mp is the Planck mass and M is the scale of infiation. It is intriguing 
to plug in the numbers. For infiation on the GUT scale one has A^pert ~ 10^ 
e-foldings before screening becomes effective. For electroweak infiation we 
predict Apert ~ 10^^ e-foldings. 

For technical and historical reasons we chose to develop the scheme for 
source 4d of Table 2: 

The dimensionless parameter e is only about 10~^^, even for GUT scale in- 
flation.0 Specializing to a homogeneous and isotropic background we have: 

F[6](t) = J^t'e-''^''^ |J'rftV*(*") (^^6(t") +36^(t")) B^[b]{t") , 

(163) 

where we define the functional B[b]{t) as follows: 

fi[6](t) =yJ^*dtV3^(*')^*'rftV(*") (^h{t") + 3b\t")^ . (164) 
The time derivative of F[b]{t) is simple to compute: 

' J^dt'e-''^'^+'''^''^ (^h{t') + 362(t')) B^[b]it') . (165) 



dF 172 ^ /•* , 



minor point is that the parameter e characterizes aU models, not just those with 
three factors of the Ricci scalar. Even when the explicit factors of A are replaced by i?'s 
one still gets when the co-moving time is rescaled to the dimensionless variable r = Ht. 
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Although the evolution equation we have obtained is both non-local and 

non-linear, it is simple to solve numerically by naive discretization. The 
independent variable t is characterized by an integer i and a dimensionless 
step size At: 

t — > iH-^Ar . (166) 
Functions of t become discrete in the usual way: 

fit) fi . (167) 

Derivatives are discretized using the difference operator: 

m H = . (168) 

In order to preserve the fundamental theorem of integral calculus we must 
discretize integrals by summing one step backwards: 



fdt'fit') 
Jo 



E^-'Ar/,-- (169) 

j=0 



It turns out that all the factors of H ^ Ar cancel for our ansatz (this is its nice 
technical feature) and we obtain the following discretized evolution equation: 

A^,.-i(^V, (170) 



2 Vl-F, 

where the discrete versions of F[b]{t) and B[b]{t) are: 
172 '-^ 



Fi ^ ^e'j:e-''^j:e''''(^A\ + 3{Ahr){B,)\ (171) 

j=0 fe=0 ^ 



j=Q k=0 

Bi = 

j=0 k=0 



We have simplified the notation by dropping the hat on the discretized version 
of -F[&](t) because F[g]{t,x) and F{t) are never discretized. 

Although the discretization we have achieved is plausible, it would be 
tedious to iterate on account of the need to compute summations at each 
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step. This can be avoided by simply keeping AFj and ABi as auxiliary 
variables. The resulting recursion scheme is: 



Bi = 5i_i + A5,_i, (173) 

AB, = 6-=^^"-^ (a5,_i + ^AVi + 3(A6,_i)2) , (174) 

= i^,_i + AF,_i, (175) 

AF, = e-=^^''-(AF,_i + lI^62gA%_i + 3(A6,_0')(5.-i)') ,(1^^^ 

k = 6.-i + A6,_i, (177) 

Ab, = A6i_i + AVi, (178) 

= -U-^)' ■ (179) 



2 Vl-F, 

Note that only differences of b enter the scheme; the exponentials of ±36 in 
the continuum theory are all absorbed into such terms. Note also that all 
the variables are initialized to zero except Afog = At. It is only through this 
initial value that the scheme depends on the step size At. 

Although e should be about 10~^^ or less, the evolution is very slow for 
such small values. In this first study we accordingly ran the scheme for 
somewhat larger, although still sub-Planckian values: e = 10~^, 10~^, 10~^ 
and 10~^. The evolution was carried out with a step size of At = 10~^ 
using the package Mathematica [^]. The results for the effective Hubble 
constant -ffefr(^) = b(t) are displayed in Figure 1. In each case the end of 
inflation is rapid and requires about five e-foldings, in good agreement with 
the perturbative prediction 0. 

Figure 2 shows the instantaneous equation of state ptot/ptot for the four 
values of e. The asymptotic equation of state is clearly that of pure radiation: 
Ptot = |ptot, corresponding to a scale factor which grows as the square root 
of the co-moving time. The reason for this is that the functional F continues 
to approach one as long as i? = 6(6 + 26^) is positive. As F approaches one, 
the second derivative of b{t) become ever more negative, which drives R to 
zero. But this implies square root expansion: 

R{t)=0 =^ 6(t) = _l_, (180) 

where the shift provides a reasonable definition for the time at which 
inflation ends. 
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Effective Hubble Constant for eps = .001 



Effective Hubble Constant for eps = .0001 
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Equation of state for eps = .001 



Equation of state for eps = .0001 
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Figure 2: ptot/ptot versus Ht for e = 10-^ 10"^, 10"^ and 10"^ 
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e 


M (GeV) 


-^pert 




a 


Htp 


10-3 


1.7 X 10^^ 


37.4 


38.1 


.2085 


.04959 


10-4 


9.6 X 10^^ 


173.6 


174.9 


.2347 


.01208 


10-^ 


5.3 X 10^^ 


805.8 


808.7 


.2425 


.002680 


10-*^ 


3.0 X 10^^ 


3740. 


3748. 


.2447 


.0005837 



Table 3: Parameters from the various runs. 



Of course the hmiting form is approached asymptotically. The first cor- 
rections in the series are: 

m = 1-^ + ... , (182) 

where the higher corrections are down by inverse powers of — t^), possibly 
offset by logarithms of same. The parameters Htz, a and Htp are determined 
by fitting the curves for each of the four runs. Their numerical values are 
reported in Table 3. Points to note are the close agreement between iVpert 
and Htz, and the near constancy of a over three decades of variation in e. 
The first fact means that the number of e-foldings for inflation to end is well 
predicted by perturbation theory; the second fact means that the transition 
to radiation domination is almost independent of the scale of inflation. That 
the transition is also quite rapid is illustrated by Figure 3. Note that it 
is extremely rapid with respect to the evolving time constant provided by 
Hes{t) because the corrections fall off with the time constant H ^ Hcs{t). 

Once square root expansion is accepted one can actually derive the asymp- 
totic series from the evolution equation. To see this, substitute the asymp- 
totic forms: 

Pmt) = l-SF(t), (184) 

into the evolution equation ( |159| ) . Taking the square root gives a differential 
equation for 6F(t): 

I ln(5F) = -^_+it- tz)m + . . . (185) 
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Pure radiation (solid) versus numerical (dashed) for eps = .000001 



Asymptotic formula (solid) versus numehca! (dashed) for eps = .000001 





Figure 3: Asymptotic approach to radition domination for e = 10~^. Left- 
hand graph shows pure radiation (sohd) versus the numerical result (dashed). 
The righthand graph show the asymptotic formula (solid) versus the numer- 
ical result (dashed). 



Integration gives the constant tp iia (|182|) . To get ( |181| ) we differentiate 
expression ( |165| ): 

_Sb^ + 86e^b + 2b')B^[b]{t) . (186) 



dt 

Now substitute ( |182| ) and the leading term in ( |181|) to obtain the following 
leading order result: 



df^ dt 2{t-tFf 



;i87) 



This must be equal to 86e^(6 + 26^)5^[6]. The time dependence must come 
from the term b{t) + 2b'^{t) since B[b]{t) is dominated by what went on during 
inflation: 



B[b]{t) = B[b]{Q + 
= HU + ... 



The asymptotic form ( |183| ) gives: 



b{t) + 2b\t) = 6b{t) + 



25b{t) 
t - 



+ 



(188) 
(189) 



(190) 
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Substituting everything in ( |18(j| ) gives a differential equation for Sb(t): 



^ + ._86.w(*M0-.g)+.... (191) 



The solution is the second term in ( |181| ) with the relation: 



Htp ~ n2e\Ht,Ya . (192) 

The independently fitted parameters of Table 3 obey ( |192| ) to a few per- 
cent, with agreement better for smaller values of e. The slight discrep- 
ancy is probably due to the approximation -B[&](t) ~ B\bc\ass\iiz) = Ht^i 
which is an overestimate. If we assume that a is almost constant then 
Htz ~ Alport ~ e~^/^ implies that Htp varies as the two thirds power of 
e. A consequence is that the post-inflation value of F[g\{t) must be very 
close to one for low scale inflation. 

It is perhaps significant that the effective Hubble constant h(t) approaches 
the expansion rate for a radiation dominated universe from below. Although 
matter is negligible during inflation, it cannot be ignored afterwards. We 
must therefore expect that matter radiation is produced during the tran- 
sition. The asymptotic expansion rate means that this matter radiation 
is progressively enhanced with respect to the purely quantum gravitational 
stress energy we have been discussing. 



7 The scalar potential 

The purpose of this section is to reconstruct the scalar potential for the 
model we have just evolved. We begin by reviewing the general technique, 
then we obtain analytic expressions during the perturbative and late time 
regimes. The section closes with an explicit numerical reconstruction over 
the full period of evolution. 

Recall from Section 2 that the scalar potential V{(j)) does not enter the 
evolution equation and is therefore not required to determine b{t). We instead 
reconstruct the potential from b{t) by imposing stress energy conservation. 
The procedure is first to get the potential as a function of time: 

^(t) = rr(^(t) + ^^"(t)-^H') . (193) 
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One then inverts the relation: 

= - 



ln(l - F[h]{t) 



(194) 



to find time as function of 0. Substituting in ( p.93| ) gives V^(0). 

This procedure can be carried out analytically during the perturbative 
regime where we can find explicit expressions for (|193|) and ( |194| ): 



V{t) 



A 



{Ht/N, 



pert/ 



SvrG ATpert 1 



{Ht/N, 



1 - 



{Ht/N, 



pert J 



pert ) 



4A^pert 1 



{Ht/N, 



pert ) 



V8^ 



In (1 - (Ht/N, 



pert ) 



,(195) 
(196) 



(Recall that A^pert is our perturbative estimate ( |161| ) for the number of e- 
foldings of infiation.) Solving for the time as a function of the scalar: 



Ht 



pert 



1-e- 



/8ttG<I> 



(197) 



gives the following relation for the potential during the perturbative regime: 



8nG N, 



pert 



(1 



STTG<t> 



/8ttG<I> 



4K 



pert 



(1 



/SnGip 



. (198) 



This expression should be valid for < < ln(A'pert)/v87rG- Note that it is 
the same for all models since it follows from the known perturbative results. 

The other regime where explicit expressions can be obtained is that of 
late times. For the model of Section 6 we have: 



A 



SttG 
1 



1 



:ln 



1 



+ 



V8^ \t-t^ 

During this period the time can be expressed in terms of the scalar as: 



(199) 
(200) 



t = tz + tpe^ 



SttG( 



+ 



(201) 
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The Effective Scaiar Potentiai for epsiion = .001 , .0001 , .00001 , and .000001 



? -0, 




1 1.5 

phi in units of M_Pianck 



Figure 4: V{(j)) versus for e = 10 ^ (leftmost), 10 ^, 10 ^ and 10 ^ (right- 
most). 

which gives the following result for the late time potential: 

2 



V|ate(0) = - 



A 



8nG 



1 



12 I Ht, 



+ . 



(202) 



This expression should be valid for > — \n{HtF) /VSttG. Using ( |192|) we 



see that this is about > \n{Npert) / V^nG , so most of the range is covered 
by the two asymptotic expressions. 

Full coverage can be obtained by computing the values of and V"(0) at 
each step in the numerical evolution, and then plotting the resulting curve. 
The discretized formulae are: 



Mr 



in 



ln(l-F,), 



(203) 



44 



M4 



Ah 



) 



2 



1 - 



Ar 



3 Ar2 • 



(204) 



The curves for each of the four runs are plotted in Figure 4. Note that they 
confirm the results of our asymptotic expansions, including the fact that the 
transition from inflation occurs at ~ In(A^pert) / V SnG. 

8 Discussion 

We have developed a compelling physical picture for a particular second or- 
der back-reaction that quantum gravity has on inflation. The first order 
effect is that long wavelength virtual gravitons can be ripped apart by the 
super luminal expansion of spacetime. This is the phenomenon of super adi- 
abatic amplification, first studied by Grishchuk [Q. This first order effect is 
not secular. As more and more graviton pairs are injected into the inflating 
universe, the growth of their energy is cancelled by the expansion of the 3- 
volume to produce a small, constant energy density (and pressure p = —p) 
of magnitude about H^. The secular effect comes at the next order from the 
gravitational interaction between the receding virtual gravitons. As graviton 
pairs are pulled apart their long range gravitational potentials fill the inter- 
vening space, and these potentials remain to add with those of new pairs 
even after the old pairs have been redshifted into insignificance. The sec- 
ond order effect is suppressed by the small dimensionless coupling constant 
GA < 10~^^, but it is cummulative. And it slows inflation because gravity is 
attractive. 

An explicit two loop computation has already confirmed that the quantum 
gravitational back-reaction slows infiation by an amount which eventually 
becomes non-perturbatively large 0. The question is, what happens next? 
We have argued in this paper that the question can be answered by inferring 
and then numerically evolving the most cosmologically significant terms in 
the effective field equations. Section 2 proved that the important part of the 
quantum gravitationally induced stress tensor is that of an effective scalar 
field (f)[g\ which is itself a non-local functional of the metric. We also showed 
that a model is completely specified by (t>[g\, since the potential Vi^cj)) does 
not enter the evolution equation, and we obtained an explicit expression (|30D 
for the scalar during the perturbative regime. 

In Sections 3 and 4 we showed that general considerations and the physics 
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of screening very largely constrain the scalar. It must have the form: 

0b] = -^7i^ln(l-F[^]) , (205) 

where the functional F[g\ consists of a retarded scalar Green's function D"^ 
acting on a source composed of two more factors of D"^ acting in some 
order on from three to five Ricci scalars. We did not determine whether 
the two inner Green's functions act in series or in parallel, nor did we fix 
how the various Ricci scalars are located with respect to them. The resulting 
combinatorics yields 73 possibilities, many of which were tabulated in Section 
5. We also identified two procedures for actually deriving the scalar from the 
result of perturbative computations of the same complexity as the one already 
done. 

In Section 6 we selected one of the candidate models: 

and numerically evolved it through the end of inflation into the regime of 
asymptotically late times. We found that inflation ends over a period of 
about five e-foldings, following which the universe asymptotically approaches 
the square root expansion of a radiation dominated universe: 



m=-'4mri(RaRaR]]], (208, 



Almost identical results were found for the series analog model: 

although we did not report them. We were able to confirm the asymptotic 
forms by deriving analytic expressions from the non-linear and non-local 
evolution equation. An interesting consequence of this analysis is that the 
approach is from below, implying that any matter radiation produced in the 
transition would be relatively enhanced by the slower redshift. 

The asymptotic analysis can be used to categorize models by the number 
of factors of R which lie immediately to the right of the outer n~^. Suppose 
there are / outer factors of R. We can expose them by taking derivatives: 

cPF -dP , ^ 
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The two internal retarded Green's functions will be effectively constant, dom- 
inated by what went on during inflation, so the time dependence must come 
from the / outer factors of R. Suppose that the asymptotic expansion rate 
is a general power law p > 0: 

b{t) ^ . (210) 

The Ricci scalar goes to: 

m = Qb{t) + 12b\t) ^ , (211) 

so the lefthand side of ( p09| ) goes like (t — tz)^"^-^, unless p = 1/2. To find the 
righthand side, consider the evolution equation: 

It follows that F[b]{t) must have the form: 

F[6](t)-^l-(^j • (213) 



Hence the righthand side of (209) goes like (t — tz) ^ and we must have 



2 + ^2^ = 2/, unless p = 1/2. 

It follows that the case of no outer Ricci scalars (/ = 0) is not consistent 
with stable evolution. What happens for these models is that F\b]{t) actually 
reaches one and b{t) goes to minus infinity. For / = 1 we get p = 0, which 
is also not consistent. However, taking account of next order terms in the 
asymptotic expansion of b gives p = 1/2. The other cases are all consistent: 
for / = 2 we get p = 2; / = 3 gives p = 8; / = 4 results in p = 18; and / = 5 
produces p = 32. Since the actual model is likely to be a linear combination 
of the tabulated candidates, and since the term with the lowest value of 
/ dominates the asymptotic behavior, we conclude that the approach to a 
radiation dominated universe is generic for models which include at least one 
/ = 1 term while avoiding any / = terms.Q 

^The f ~ terms on Table 1 are 16a-21a and 10b-12b; the f — I terms are lla-15a 
and 7b-9b. On Table 2 the / = terms are 7c-10c and 5d-6d; the f — 1 terms are 4c-6c 
and 3d-4d. 
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Of course fixing the asymptotic time dependence does not provide a phys- 
ical interpretation for what happens after the end of infiation. For example, 
the stress energy of our model is not likely to consist of gravitons, in spite of 
the fact that it approaches the equation of state of radiation. There is sim- 
ply no way to have created them. The energy density of gravitons produced 
by superadiabatic amplification is a small constant for as long as it can be 
reliably tracked using perturbation theory — which is almost to the end of 
infiation. What terminates infiation is the buildup of gravitational interac- 
tion stress energy and it must be this which adds with the stress tensor of 
the cosmological constant to produce a residual obeying the equation of state 
of radiation. 

The effective scalar (f)[g\ poses a similar interpretational dilemma. It is 
certainly not a fundamental particle but one might plausibly interpret screen- 
ing in terms of the formation of a scalar bound state on cosmological scales. 
The "scalars" would be the virtual gravitons pairs which are ripped apart by 
infiation. Although they are separated to cosmological distances it is their 
binding energy which eventually arrests infiation. Since the pairs tend to an- 
nihilate on sub-horizon scales there are no new massless quanta to embarrass 
phenomenology. One very attractive feature of this interpretation is that 
we can compute the spectrum of density perturbations using the standard 
formalism, just as if the scalar was fundamental. 

Finally, there is the question of how to couple matter and what effect 
doing so will have, both on the geometry and on reheating. Gravitons dom- 
inate screening through their unique combination of masslessness without 
conformal invariance, however, ordinary matter becomes important at the 
end of infiation. The simplest assumption would that the matter and the 
gravitational stress tensors are separately conserved except for a brief period 
at the end of infiation when some of the gravitational stress energy excites 
matter degrees of freedom. It remains to see if there can be sufficient reheat- 
ing without the phase of coherent oscillations that characterizes scalar based 
infiation p|, 
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